function [tempM, valueF0, vF1, difFF]=TissueS_miRNA_biomass(expressionData_tissue,miRNA_Gene, model00)

[NumMiRNA,NumTar]=size(miRNA_Gene);

[NumLocus]=size(expressionData_tissue.Locus,1);

miRNA_metaGene=zeros(NumMiRNA,NumLocus);

for i_mirna=1:NumMiRNA
    for j_tar=1:NumTar
        if isnan(miRNA_Gene(i_mirna,j_tar))
        else
        
            ind=find(expressionData_tissue.Locus==miRNA_Gene(i_mirna,j_tar),1);
            if ind
                 miRNA_metaGene(i_mirna,ind)=1;
            end
        end
    end
end

tempM=miRNA_metaGene;

[modelF0,temprxn] = createTissueSpecificModel(model00,expressionData_tissue);

modelF0=changeObjective(modelF0,'bioMassR');
valueF0=optimizeCbModel(modelF0);


for imodel=1:NumMiRNA
%for imodel=77:77    
    EXPtemp=and(expressionData_tissue.Data, not(tempM(imodel,:)'));
    difFF(imodel,1)=sum(expressionData_tissue.Data)-sum(EXPtemp);
    EXPi=expressionData_tissue;
    EXPi.Data=double(EXPtemp);
    
%    [temptissueModel,temprxn] = createTissueSpecificModel(modelF0,EXPi);
    [temptissueModel,temprxn] = createTissueSpecificModel(model00,EXPi);
    Fi=optimizeCbModel(temptissueModel);

    vF1(imodel,1)=valueF0.f-Fi.f;
    
end
return;